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Abstract 

We investigate the nonperturbative renormalisation of composite operators in lattice 
QCD restricting ourselves to operators that are bilinear in the quark fields. These 
include operators which are relevant to the calculation of moments of hadronic 
structure functions. The computations are based on Monte Carlo simulations using 
quenched Wilson fermions. 
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1 Introduction 



Monte Carlo simulations of lattice QCD have evolved from spectrum calcula- 
tions to more detailed investigations of hadron structure. These more advanced 
studies require to calculate hadronic matrix elements of composite operators. 
In order to compute the moments of hadronic structure functions, for exam- 
ple, one needs the matrix elements of composite operators appearing in the 
operator product expansion of the appropriate currents [1-4]. In general, one 
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has to renormalise these operators in order to obtain finite answers in the con- 
tinuum limit. Furthermore, comparison with the results of phenomenological 
analyses usually requires the matrix elements to be given in one of the popular 
continuum renormalisation schemes, e.g. the MS scheme. So one has to think 
about the conversion of the bare lattice operators to renormalised continuum 
operators. 

Consider the special case of the operators determining the moments of hadron- 
ic structure functions. Here the renormalised continuum matrix element has 
to be multiplied by the corresponding Wilson coefficient to yield the desired 
moment, which can be measured in deep inelastic scattering experiments. Be- 
ing an observable quantity it must be independent of any choices made in 
the renormalisation procedure. This is of course only possible, if the renor- 
malised matrix element and the Wilson coefficient are calculated in the same 
scheme. In particular, the dependence on the renormalisation scale /x has to 
cancel between the renormalisation constant and the Wilson coefficient. Since 
the Wilson coefficients are usually computed in the MS scheme, we have to 
convert our lattice operators to continuum MS operators if we want to make 
contact with phenomenology. 

One obvious possibility to calculate the necessary renormalisation factors is 
lattice perturbation theory. However, quite often lattice perturbation theory 
seems to converge rather slowly. Identifying one source of these poor conver- 
gence properties, Lepage and Mackenzie proposed as a remedy the so-called 
tadpole improved perturbation theory [5]. Still, lattice perturbation series 
rarely extend beyond the one-loop level, and hence considerable uncertainty 
remains. 

It is therefore only natural to try a nonperturbative renormalisation by means 
of Monte Carlo simulations. A way how to do this was introduced in Ref. [6]. 
We shall present several improvements of this method and apply it to a va- 
riety of operators that are bilinear in the quark fields [7]. In particular, we 
study operators which are needed to calculate hadronic structure functions. 
Computing the Z factors for a rather large range of renormalisation scales we 
can find out at which scales (if at all) perturbative behaviour sets in such that 
the multiplication with the perturbative Wilson coefficients makes sense. 

In this paper we shall work only with Wilson fermions in the quenched ap- 
proximation. However, an obvious next step is the use of improved fermions in 
order to reduce cut-off effects. Like the renormalisation of our operators, the 
improvement should be nonperturbative. As far as the action is concerned, 
there is already a lot of experience how to achieve this. But improvement of 
the action is not sufficient, the operators have to be improved as well. Here a 
good deal of work still has to be done, in particular for operators containing 
derivatives. 
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The paper is organised as follows: After introducing in Section 2 the opera- 
tors to be studied we explain the method of nonperturbative renormalisation 
in Section 3. In the special case of the vector and axial vector currents we 
prefer a variant of this method because it seems to suppress lattice artifacts 
more efficiently. We describe it in Section 4. Perturbative formulae, which we 
need for comparison and for obtaining results in the MS scheme, are given in 
Section 5. Section 6 contains some details of the numerical implementation 
including our momentum sources which greatly reduce the statistical noise of 
the data. In Section 7 we present and discuss our results. Finally, we come to 
our conclusions. 



2 The operators 

In the Euclidean continuum we should study the operators 

{ «l» n = qi^D, 2 ... D, n q, (2.1) 

or rather 0(4) irreducible multiplets with definite C-parity. In particular, we 
obtain twist-2 operators by symmetrising the indices and subtracting the 
traces. In the flavour-nonsinglet case they do not mix and are hence mul- 
tiplicatively renormalisable. 

Working with Wilson fermions it is straightforward to write down lattice ver- 
sions of the above operators. One simply replaces the continuum covariant 
derivative by its lattice analogue. However, 0(4) being restricted to its finite 
subgroup H(4) (the hypercubic group) on the lattice, the constraints imposed 
by space-time symmetry are less stringent than in the continuum and the pos- 
sibilities for mixing increase [2,3,8,9]. Guided by the H(4) classification given 
in Ref. [8] we have chosen the operators [4,10,11] 
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They are labeled by the reduced hadron matrix elements which they deter- 
mine. In addition, we have studied the following operators without derivatives 
("currents"): 



S = qq, (2.12) 

O p = q l5 q, (2.13) 

C£ = ffVZ, (2-14) 

(^ = ^75?. ( 2 - 15 ) 

where all quark fields are taken at the same lattice point, as well as the con- 
served vector current 

Jp(x) = \(q{x + /i)(7 M + l)U + (x, fi)q(x) 

+ g(x)( 7M - l)U(x, fi)q(x + /})) (2.16) 

with the link matrix U(x,[x) G SU(3) representing the gauge field. 

Note that the operators V2a and V2b although belonging to the same ir- 
reducible 0(4) multiplet transform according to inequivalent representations 
of H(4). Hence their renormalisation factors calculated on the lattice have no 
reason to coincide. The same remark applies to O r2a and O r2b . The operators 
O ai and O r . 2a) on the other hand, are members of the same irreducible H(4) 
multiplet, hence their renormalisation factors should agree also on the lattice. 

Concerning the mixing properties a few remarks are in order. Mixing with 
operators of equal or lower dimension is excluded for the operators V2 , V2b , 
O ai) O a2) O r2 a) O r , 2b , as well as for the currents. The case of the operator V3 , 
for which there are two further operators with the same dimension and the 
same transformation behaviour, is discussed in Refs. [8,9]. The operators V4 , 
O r:i , on the other hand, could in principle mix not only with operators of the 
same dimension but also with an operator of one dimension less and different 
chiral properties. It is of the type 

<^-> <^-> <^-> 
9<V75 D k D v*'~ D^q , ( 2 - 17 ) 



where n = 2 in the case of V4 and n — 1 for O r3 . 
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Our analysis ignores mixing completely. This seems to be justified for V3 . 
Here a perturbative calculation gives a rather small mixing coefficient for one 
of the mixing operators [3,9], whereas the other candidate for mixing does not 
appear at all in a one-loop calculation of quark matrix elements at momentum 
transfer zero, because its Born term vanishes in forward direction. The same is 
true for all operators of dimension less or equal to 6 which transform identically 
to Vi : Their Born term vanishes in forward matrix elements, hence they do 
not show up in a one-loop calculation at vanishing momentum transfer. In 
the case of O r3 , however, the mixing with an operator of lower dimension is 
already visible at the one-loop level even in forward direction. 



3 The method 



Our calculation of renormalisation constants follows closely the procedure pro- 
posed by Martinelli et al. [6]. It mimics the definitions used in (continuum) 
perturbation theory. We work on a lattice of spacing a and volume V in Eu- 
clidean space. For a fixed gauge let 

Gap{p) = ^ E e- w ^W(#)> (3.1) 

x,y,z 

denote the non-amputated quark-quark Green function with one insertion of 
the operator O at momentum zero. It is to be considered as a matrix in colour 
and Dirac space. With the quark propagator 

Sapip) = £Ee-*(^><<z«(*)fc(y)> (3.2) 

x,y 

the corresponding vertex function (or amputated Green function) is given by 

T(p) = S- 1 (p)G(p)S- 1 (p). (3.3) 

Defining the renormalised vertex function by 

r R (p) = Z q l Z Y{p) (3.4) 

we fix the renormalisation constant Zq by imposing the renormalisation con- 
dition 

^tr (rR^rBorn^)- 1 ) = 1 (3.5) 
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at p 2 = fi 2 , where \i is the renormalisation scale. So we calculate Zq from 

z- x z fcx (r(p)rBorn(p)- 1 ) = i (3.6) 



with p 2 = fi 2 . Here rB or n(p) is the Born term in the vertex function of O com- 
puted on the lattice, and Z q denotes the quark field renormalisation constant. 
The latter can be taken as 

qKF) 12Exsin 2 (ap x ) 1 } 



again at p 2 = /x 2 . 

If the operator under study belongs to an 0(4) multiplet of dimension greater 
than 1, i.e. if it carries at least one space-time index, the trace in Eq. (3.6) 
will in general depend on the direction of p. So the renormalisation condition 
(3.5) violates 0(4) covariance even in the continuum limit. In the continuum, 
this disease is easily cured by a suitable summation over the members of the 
0(4) multiplet. On the lattice, this makes sense only in the few special cases 
where the 0(4) multiplet is irreducible also under the hypercubic group H(4), 
because the renormalisation constants for the different H(4) multiplets within 
a given 0(4) multiplet will in general be different. So we have to live with 
this noncovariance, which, of course, should finally be compensated when we 
convert our results to a covariant renormalisation scheme. 

The scale \i at which our renormalisation constants are defined should ideally 
satisfy the conditions 

1/L 2 « A 2 « /i 2 « 1/a 2 (3.8) 



on a lattice with linear extent L. Whether in a concrete calculation these 
conditions may be considered as fulfilled remains to be seen. 



4 A special case: vector and axial vector currents 



In the special case of the vector and axial vector currents in the continuum we 
can distinguish between longitudinal and transverse components with respect 
to the momentum. One may thus define two renormalisation constants, one 
from the longitudinal and one from the transverse components. Denoting the 
vertex function of the currents generically by we have for the vector 
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current (the modifications required in the case of the axial vector current are 
obvious) 



Z-'Z^tr ( £ ( 7 „ - p A J,{ P )\ = 3 . (4.2) 



Note that these renormalisation conditions are 0(4) covariant. 

We try to generalise (4.2) to the lattice imposing a renormalisation condition 
of the form 

Z^Zj^tr (£ QiUp)^ = £tr QtJ B ™(P)) = 3 » ( 4 - 3 ) 



where Qj^ is a suitable Dirac matrix and J^ OTn {p) denotes the Born term 
evaluated on the lattice. In the case of the conserved vector current we know 
that Zj — 1, hence we can calculate Z q from (4.3). This value is then used to 
determine Zj for the local vector and axial vector currents. 

The matrix Qj[ is constructed as 

^ = (^t 1 - ^ > ( 4 - 4 ) 

where 

Ql = 1 sin(ap M ) J2 (^ B ° rn ) _1 Mapx) ■ (4.5) 



The normalisation factor N follows immediately from the condition (4.3): 



N = ±tr 



E ^ B ° m ^P,) E (^A B ° rn ) _1 sin(ap A ) 



(4.6) 



Note that for a — > and J^ orn — > 7 M we have iV/a 2 — > p 2 , and (4.3) reduces 
to (4.2) in the continuum limit. 

In the case of the conserved vector current we find 

J* orn (p) = 1, cos(ajv) + i sin(a^) (4.7) 



7 



so that 



N = [Y sin2 ( a p^) + Y s[n2 ( a Pv) - Y sin4 ( a Pn) ■ 



(4.8) 




For the local vector and axial vector current we have 



J, 



rBorn 



= 7m 



(4.9) 



and 



rBorn 



(4.10) 



respectively, and we obtain in both cases 



N = ^sin 2 (ap M ) . 



(4.11) 



Working with the longitudinal component, i.e. with a lattice version of (4.1), 
leads to a less smooth momentum dependence, i.e. to stronger lattice effects. 



5 Input from perturbative calculations 

Eq. (3.5) defines a renormalisation scheme of the momentum subtraction type, 
which we call MOM scheme. Note that it will in general not agree with any of 
the momentum subtraction schemes used in continuum perturbation theory. 
It is desirable to convert our results into a more popular scheme like the MS 
scheme. Moreover, we want to use our renormalisation factors in connection 
with the Wilson coefficients, which appear in the operator product expansion, 
and these are generally given in the MS scheme. Hence we have to perform 
a finite renormalisation leading us from the renormalisation scheme defined 
by Eq. (3.5) to the MS scheme. The corresponding renormalisation constant 
-^mom i s computed in continuum perturbation theory using dimensional reg- 
ularisation. 

We work in a general covariant gauge (gauge parameter £) such that the gluon 
propagator has the form 




(5.1) 
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The Landau gauge, which is employed in our numerical simulations, corre- 
sponds to £ = 0. Most of our operators can be written in the form 



J2 c mi-Mu ( ?7{m Di*"' D ^}lT 0( l ~ tracesj (5.2) 



where n 5 — 0, 1 and c^...^ is totally symmetric and traceless. Neglecting 
quark masses one obtains with an anticommuting 75 



„2 

7 ms _ 1 , y_n 

^mom — 1 + — o W 



16tt 2 



G„ + (1-0^_i 



+ 



(-^T+a-f) 2 )^ 

V n + 1 nj p 2 



2 



p 2 E M Mp) 

where CV = 4/3 for the gauge group SU(3), 



2 



+ <V), (5.3) 



2 

Cn — -7 — : 7T ( — 3 — S n ^i + 2S n+ i) 
n(n + 1) 

2 n 1 

+ r-T- 4 E-(^-Vi)-i. (5.4) 

' j=2 J 

n 2 

^n = E-> ( 5 - 5 ) 
J'=l 

and 

fyi(P)= E C W2...MnPM2---^n- (5-6) 



Because of the noncovariance of our renormalisation condition, Z^q M depends 
on the direction of the momentum p and on the coefficients c Ml) ... j/ln . 

For and C p we obtain 

= 1 + ^(4 + OCf + 0(g*) . (5.7) 

In this special case we can go one step further and use the two-loop result for 
^mom [12]. For three colours and nf flavours one has in Landau gauge: 

ZSm = 1 + f £- 2 + (177.48452 - £ „,) + <V) . (5.8, 
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Table 1 

Finite contributions to the renormalisation factors in lattice perturbation theory. In 
addition, the one- and two-loop coefficients of the anomalous dimension are given 
for flavours. 
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A 

Ao 


7o 


71 




1.2796 


64/9 


96.69 




2.5619 


64/9 


96.69 




-12.1274 


100/9 


141.78 




-27.2296 


628/45 
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-12.1171 
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o n , a 


0.3451 


64/9 


96.69 


Or 2 , b 


0.1674 


64/9 


96.69 


Or, 


-12.8589 


100/9 


141.78 


Q S 


12.9524 


-8 


-134.67 


Q P 


22.5954 


-8 


-134.67 


qV 


20.6178 





0.00 


qA 
M 


15.7963 





0.00 



We want to compare our nonperturbative results with the corresponding values 
obtained in (tadpole improved) perturbation theory on the lattice. For the 
renormalisation factor which brings us from the bare lattice operator to the 
renormalised operator in the MS scheme, lattice perturbation theory yields 
results of the form 

Z o n = 1 " 77^(70 InM + C F A ) (5.9) 



where is a finite constant and 70 is the one-loop coefficient of the anomalous 
dimension. Working with an anticommuting 75 also in the continuum part of 
the calculation we arrive at the values given in Table 1 (see Ref. [9] for more 
details). 

In order to obtain the corresponding results in tadpole improved perturba- 
tion theory [5] we write (with /i = 1/a) for an operator with n D covariant 
derivatives 
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where 



"o 



9_ 

16tt 2 ' 



-,*2 



16tt 2 



C F A + 0(g") , 



(5.10) 



«o = (|trC/ D )5 = 1 - y^CVtt 2 + 0(^ 4 ) 



(5.11) 



and 



A = A + (ra D - 1)tt 2 



(5.12) 



This reflects the fact that one has no operator tadpole diagrams and one 
leg tadpole diagram, which are of the same magnitude but contribute with 
opposite sign. Furthermore, one chooses as the expansion parameter g* the 
coupling constant renormalised at some physical scale. We have taken g* from 
the values given for a^g(l/a) in Table I of Ref. [5]. 

At this point we have two options. Either we stay with the expression (5.9) 
and its tadpole improved analogue 



Lin 



o ~ u 



1-71D 



1 



9 



*2 



16tt 2 



(7 ln(a/i) + C F A ) 



(5.13) 



or we apply these formulae only at a fixed scale jj, = fi (e.g. fi = 1/a) using 
the renormalisation group to change \i. 

For this scale dependence, the renormalisation group predicts (at fixed bare 
parameters) 



ROW, 1*0) ■= „ 7 x = eX P \ 



9(M 2 ) 



dg 



jig) 



> . 



(5.14) 



Here the /5-function is given by 



P(g) = -A) 



16tt 2 



Pi 



'16tt 



2\2 



+ 



(5.15) 



with 



2 38 
ft = 11 -- n/ , ft = 102 



(5.16) 
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In terms of the A parameter the running coupling reads 

flV) = 1 _ Pi mln(/x 2 /A 2 ) 

16tt 2 p ln(/i 2 /A 2 ) PI ln 2 (/i 2 /A 2 ) 

and the anomalous dimension is expanded as 



+ ••• (5.17) 



^HTo^+J^V-- (5.18) 



For our operators the one- and two-loop coefficients 70 and 71 of the anomalous 
dimension are known in the MS scheme [13,14]. The values actually used are 
listed in Table 1. Within the two- loop approximation we obtain 



/l j_ /3 1§ 2 (M 2 )\ ^ ({7l//3l) " {7 ° //30)) 



1 + ft 16tt 2 

1 , ft g 2 (^g) 
\ 1 ^ ft 16tt2 



(5.19) 



In the case of the operators O s and O p we can do better and use the three- 
loop expressions for the P function and the anomalous dimension. In the MS 
scheme one finds [15-17] 



n 2857 5033 325 2 

p 2 = rif H n f , 

f 2 18 ; 54 J ' 

404 4 

7i = -— + (5.20) 
/4432 320 7 280 2 



Z4432 320 A \ 
72 = -2498 + l^— + — C 3 J % + 



81 f 



where (3 ~ 1.2020569032. For the running coupling we have now instead of 
(5.17) " 



tfV) 1 Alnln(/x 2 /A 2 ) 



16tt 2 A,hi(/i 2 /A 2 ) pi InV/A 2 ) 

+ IhPo - ft 2 ) + • • ' (5.21) 

and i?o takes the form 

Ro{p,IM>) = y&y{ (5-22) 
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with 



F(x) = £ Inx + In (Po + Pi* + fk* 

+ VMto-frbrH-MT, arctan ( ft + 2ft* , (5 23) 



Having data at two different values of the bare coupling we can also compare 
the dependence on the bare coupling with the behaviour expected from the 
renormalisation group. Keeping the renormalisation scale /x and the renor- 
malised quantities fixed we find for the ratio of the renormalisation factors Zq 
and Z' at the bare couplings g and g', respectively, 



Z'q 



o 

T = exp < 



9 A 7(0o) 
dg - 



(5.24) 



Here 7(g) and /3(g) denote the anomalous dimension and the /^-function ob- 
tained by differentiation with respect to the cut-off at fixed renormalised quan- 
tities. They are considered as functions of the bare coupling. Within the two- 
loop approximation we get 




-l((7l/A)-(70/A))) 



Z' VJ -j \ ! 



(5.25) 



where we have used the fact that p , Pi, and 70 are universal so that /3 = (3 , 
Pi = Pi, 7o = To- The coefficient 71 reads (cf. also [18]) 

7l = 7l + 2p C F A o - 327T 2 7o (-0.234101 + n f ■ 0.0034435) (5.26) 



for the gauge group SU(3). The quantities Ao are given in Table 1, and the 
coefficient of 7 o stems from the ratio of the A-parameters on the lattice and 
in the MS scheme (see, e.g., [19]). 



6 Numerical implementation 



Let us sketch the main ingredients of our calculational procedure. To simplify 
the notation we set the lattice spacing a = 1 in this section. In a first step gauge 
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field configurations are generated and numerically fixed to some convenient 
gauge (the Landau gauge in our case). The non-amputated Green function 
(3.1) is calculated as the gauge field average of G(p), which is constructed from 
the quark propagator S(p) on the same gauge field configuration according to 

G(p) = ^ E e-*l*-y)S(x,z)J(z,2>)S(2>,y). (6.1) 

x,y,z,z' 

We omit the quark-line disconnected contribution (which would contribute 
only to flavour-singlet operators) and suppress Dirac as well as colour indices. 
The operator under study is represented by J(z, z')\ 

520(z) = 52q(z)J(z,J)qW. (6.2) 

z z,z' 

Using the relation 

S(x,y)= l5 S(y,x) + l5 (6.3) 

we rewrite G{p) as 

G{p) = E 75 f E S(z, x)#>°) l5 J(z, z') ( E S(z', y)&A (6.4) 

z,z> \ x ) \ y ) 

in terms of the quantities 

J2S(z,x)e^. (6.5) 

X 

These can be calculated by solving the lattice Dirac equation with a momen- 
tum source: 

E M(y, z) (E S(z, x)e^A = e^ . (6.6) 

2 \ X J 

Here M(x, y) represents the fermion matrix. So the number of matrix inver- 
sions to be performed is proportional to the number of momenta considered. 
But the quark propagators, which we need for the amputation and the compu- 
tation of the quark wave function renormalisation, are immediately obtained 
from the quantities already calculated. 

Another computational strategy would be to choose a particular location for 
the operator, i.e. instead of summing over z' (and z) in (6.1) one could set 
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z' = 0. Translational invariance tells us that this will give the same expectation 
value after averaging over all gauge field configurations. For this method we 
need to solve the Dirac equation with a point source at the location of the 
operator and (in the case of extended operators) for a small number of point 
sources in the immediate neighbourhood. For operators with a small number of 
derivatives the point source method would require fewer inversions, but it turns 
out that relying on translational invariance leads to much larger statistical 
errors. This is clearly seen in the results shown in Ref. [20], which were obtained 
with a point source instead of a momentum source (cf. also Ref. [21]). 

We work with standard Wilson fermions (r = 1) in the quenched approxi- 
mation. At /3 — 6.0 we have analysed 20 configurations on a 16 3 x 32 lattice 
for three values of the hopping parameter, k = 0.155, 0.153, and 0.1515. 
Four further configurations at f3 — 6.2 on a 24 3 x 48 lattice were studied 
with k = 0.152, 0.1507, and 0.1489. For the critical k we take the values 
k c = 0.157211 (p = 6.0) and k c = 0.153374 {f3 = 6.2), respectively (see 
Ref. [22]). 

Before we calculate quark correlation functions on our configurations we fix 
the gauge to the Landau gauge [23]. The gauge fixing necessarily raises the 
question of the influence of Gribov copies. Fortunately, an investigation of 
this problem in a similar setting indicates that the fluctuations induced by 
the Gribov copies are not overwhelmingly large and may be less important 
than the ordinary statistical fluctuations [24] . Therefore we decided to neglect 
the Gribov problem for the time being, but a careful study in the case at hand 
is certainly desirable. 

Coming finally to the choice of the momenta, we tried to avoid momenta along 
the coordinate axes in order to minimise cut-off effects. Altogether we used 43 
momenta at f3 — 6.0 and 41 momenta at f3 — 6.2. 



7 Results 



For the presentation of our results we convert lattice units to physical units 
using a~ 2 = 3.8 GeV 2 ((3 = 6.0) and a~ 2 = 7.0 GeV 2 {(5 = 6.2) as determined 
from the string tension [22]. This leads to values for the spatial extent of our 
lattices of 1.6 fm (f3 = 6.0) and 1.8 fm (f3 = 6.2). The A parameter is taken 
to be A^g = 230 MeV. This value follows from the result given in Ref. [25] 
by using the string tension instead of the force parameter r to set the scale. 
As we work in the quenched approximation we have to put rif = in the 
perturbative formulae. The plotted errors are purely statistical and have been 
calculated by a jackknife procedure. 
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Fig. 1. Chiral extrapolation of Z for V2i 
of /j, 2 have been selected. 



at = 6.2. Some representative values 



The quark masses used in our simulations are rather large. On the other hand, 
the perturbative calculations, which we use for comparison and for converting 
to the MS scheme, are performed at vanishing quark mass. Hence we first 
extrapolate our results to the chiral limit. This is done linearly in In 
most cases the mass dependence is quite weak so that the extrapolation looks 
reliable, see Fig. 1 for an example (operator V2 a at (5 — 6.2). This is not 
too surprising since the mass scale in the calculations is not set by the quark 
mass, but by the (off-shell) momentum, which is typically much larger. The 
only exception is the pseudoscalar density where a rather strong mass de- 
pendence is found. This observation already indicates that the pseudoscalar 
density is something special, and we will discuss this operator in more detail 
in Subsection 7.2. 



In Figs. 2 - 4 we display our "raw" results, i.e. the nonperturbative Z's in our 
MOM scheme, for f3 = 6.2. Then we multiply by ^mom i n order to obtain 
the Z factors in the MS scheme. When evaluating the perturbative expression 
for ^mom we insert for the coupling constant the running coupling in the MS 
scheme taken at the renormalisation scale fx 2 (cf. Eq. (5.17)). Fig. 5 illustrates 
the effect of the conversion to the MS scheme for the operator V2 a at (3 = 6.2. 

Ultimately we are interested in physical observables, e.g. in the moments of 
the structure functions. They are obtained from the bare matrix elements cal- 
culated on the lattice after multiplication by the appropriate Z factor and the 
Wilson coefficient (both computed in the MS scheme). So the bridge from the 
lattice to continuum physics consists of the product of these two quantities. 
Having calculated Z nonperturbatively we can multiply with the correspond- 
ing Wilson coefficient in order to see if the expected cancellation of the \i 
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Fig. 2. Z in the MOM scheme for <D%, Ofi, O p , and O s (from top to bottom) 
at (3 = 6.2. For the local vector current and the axial vector current, Z has been 
determined from the transverse components (cf. Section 4). 
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Fig. 5. Z for the operator V2a at /3 = 6.2. The filled circles represent the data cor- 
responding to the MOM renormalisation condition (3.5), the open circles represent 
the values in the MS scheme. 

dependence really occurs. If it does, the \i independent value of the product is 
the factor we need in order to calculate the moments of the structure functions 
from the lattice results. Since this cancellation means that the scale depen- 
dence of our Z's should be described by the renormalisation group factor R 
(given by (5.19) in two-loop approximation) we shall divide our numerical re- 
sults by this expression and define Zrqi = Z/R. For Zrgi we hope to obtain 
a n independent answer, at least in a reasonable window of ji values. We shall 
choose $ = 4 GeV 2 in (5.19). 



7.1 Vector and axial vector currents 

Let us begin with the vector and axial vector currents. Since in these cases 
the MS anomalous dimensions vanish, the renormalisation group factor (5.19) 
equals 1 and we expect to find immediately scale independent results in a suit- 
able window. As already mentioned, the standard procedure from Section 3 
suffers more strongly from lattice artifacts than the alternative method de- 
scribed in Section 4. This is exemplified in Fig. 6. The results obtained by 
applying the standard procedure to (77475(7 scatter more strongly than the Z's 
from the transverse component of the axial vector current. In the remainder 
of this paper we shall only use the values delivered by the latter procedure. 

Fig. 6 nicely fulfills our expectations. We see a flat region around /i 2 ~ 4 
GeV 2 where Z becomes scale independent as it should. Moreover, tadpole 
improvement has really improved the agreement between our nonperturbative 
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Fig. 6. Z for the axial vector current (transverse component) at = 6.0 determined 
by the method of Section 4 (upper plot). The lower plot shows Z for the operator 
a t the same (3 value. The dashed (solid) line represents the prediction of 
one-loop lattice perturbation theory with (without) tadpole improvement. 

result and one-loop lattice perturbation theory. 

The analogous plot for the transverse component of the local vector current 
is displayed in Fig. 7. No indication of a scaling window is seen. How can this 
different behaviour of local vector and axial vector current be explained? A 
hint may be obtained from lattice perturbation theory. If one keeps all terms of 
order am one gets an additional contribution of the qualitatively correct form 
in the case of the local vector current, whereas the analogous contribution to 
the renormalisation of the axial vector current is considerably smaller [26]. 
Due to the spontaneous breakdown of chiral symmetry one should expect 
that some mass-like effects survive even in the chiral limit and may lead to 
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Fig. 7. Z for the local vector current (transverse component) at (3 = 6.2 determined 
by the method of Section 4. The dashed (solid) line represents the prediction of 
one-loop lattice perturbation theory with (without) tadpole improvement. 
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Fig. 8. Z for the local vector current (transverse component). The open circles 
(filled squares) represent the data for j3 = 6.2 ((3 = 6.0). The (3 = 6.2 data have 
been rescaled perturbatively to (3 = 6.0 by multiplication with (5.25). 

the observed behaviour of the local vector current through such am terms. 

This interpretation is supported by a comparison of the data at /3 — 6.0 
and (3 = 6.2. The perturbative renormalisation group tells us that the (3 
dependence of the Z's is given by the /z 2 independent factor (5.25) in terms of 
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Fig. 9. Z for the axial vector current (transverse component). The open circles 
(filled squares) represent the data for (3 = 6.2 ((3 = 6.0). The (3 = 6.2 data have 
been rescaled perturbatively to [3 = 6.0 by multiplication with (5.25). 

the bare coupling constants g and g'. Therefore we multiply the f3 = 6.2 data 
with the factor (5.25) and plot them together with the Z's at f3 — 6.0 versus 
the renormalisation scale fj, 2 . The results are shown in Fig. 8 for the local vector 
current and in Fig. 9 for the axial vector current. In both cases, the f3 = 6.2 
data are closer to our expectations, but the local vector current seems to be 
more sensitive to the variation of a than the axial current, in accord with the 
above mentioned stronger influence of am corrections in lattice perturbation 
theory. Indeed, comparing the slopes of the data for the local vector current 
at the two f3 values in the region below fi 2 ~ 15 GeV 2 one can estimate a ratio 
of about 1.45, close to the ratio a((3 = 6.0)/a(/3 = 6.2) « 1.36. 



7.2 Scalar and pseudoscalar density 

Fig. 10 shows our results for the Z factor of the scalar density. The upper part 
of the figure compares the nonperturbative numbers transformed to the MS 
scheme with the predictions from lattice perturbation theory with and without 
tadpole improvement. We also display the curves resulting from renormalisa- 
tion group improvement applied to the perturbative results at /i 2 , = 4 GeV 2 . 
Tadpole improvement works in the right direction and the data seem to follow 
the trend required by the renormalisation group if fi 2 is not too small. This is 
more clearly exhibited in the lower part of the figure where we have divided 
by the renormalisation group factor (5.19). As expected, we find a nice scaling 
window in Zrgi although some lattice artifacts are clearly visible. 
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Fig. 10. Z and ^rgi for the scalar density at f3 = 6.2. The dashed (solid) straight 
line represents the prediction of one-loop lattice perturbation theory with (without) 
tadpole improvement (Eqs. (5.13), (5.9), respectively). The dotted curve results 
from improving the perturbative prediction with the renormalisation group (cf. Eq. 
(5.14)) and matching at $ = 4 GeV 2 . The dash-dotted curve represents the same 
modification for tadpole improved perturbation theory. 



The results for the pseudoscalar density do not agree with the naive expecta- 
tion that Z RG i should be constant (see Fig. 11). This observation was already 
made by Martinelli et al.[6]. Moreover the data show a strong mass depen- 
dence, especially for the lower values of fi 2 . These features can be explained by 
the dynamics of chiral symmetry breaking [6] . Indeed, in the continuum there 
is a Ward identity relating the (bare) pseudoscalar vertex T p (p) at momentum 
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Fig. 11. Zrgi for the pseudoscalar density at (3 = 6.2. 
transfer zero to the quark propagator: 

- 2m T p (p) = S-\ P )k + l 5 S~\p) . 



(7.1) 



Here m is the bare quark mass, which is taken to be nonzero. According 
to Pagels [27], spontaneous breakdown of chiral symmetry produces a mass- 
like contribution to S* -1 which commutes (rather than anticommutes) with 
75 and persists in the chiral limit. Thus the r.h.s. of (7.1) will not vanish 
for m — > and the pseudoscalar vertex must diverge. Consequently, the 
pseudoscalar density is ill defined if the chiral limit is performed after the 
momentum transfer has been set to zero, which is the order of limits in our 
procedure. 



For our Z factors the identity (7.1) would lead to 



Z p(p) = 



12m Z q (p) 
tr S-^p) 



(7.2) 



and with tr S~ 1 {p) remaining nonzero even in the chiral limit Zqp(jj) should 
vanish for m — > 0. This tendency is indeed shown by our data (at least for 
the lower values of ji 2 ). 

Nevertheless, for sufficiently large /x 2 , i.e. in the perturbative regime, the ratio 
Zqp/Zqs should become constant, because O p and O s have equal anomalous 
dimensions. Comparison of the lower part of Fig. 10 with Fig. 11 shows that 
our data do not display this behaviour. The ratio continues to rise up to our 
largest values of fj 2 . 
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Fig. 12. Z RG1 for O v%a at (3 = 6.2. 
7.3 Operators with derivatives 



In this subsection we want to discuss some selected operators with derivatives, 
which determine moments of hadronic structure functions. Again, we would 
like to find a "window" at moderate values of fi 2 where -Zrgi is constant. 
There /z 2 would be large enough to allow for perturbative scaling behaviour 
and the multiplication with perturbative Wilson coefficients would make sense. 
On the other hand, /i 2 should be small enough to avoid strong cut-off effects. 
Indeed, as we have seen, the results for the scalar density are rather close to 
this ideal scenario. For operators with derivatives the situation is however less 
favourable. 



Looking at the results for Z RG i one observes that one gets similar values for all 
operators containing the same number rip of derivatives. They increase with 
we. Qualitatively, this behaviour is reproduced by lattice perturbation the- 
ory, and tadpole improvement enhances the effect, though not in quantitative 
agreement with the nonperturbative data. In order to facilitate the compar- 
ison of the different operators we have chosen the same vertical scale in the 
plots for all operators with the same number of derivatives although in some 
cases this has the consequence that the results for a few low values of /i 2 do 
not appear on the plot. 

In Figs. 12-15 we plot Z RGI versus the renormalisation scale fx 2 for various 
operators containing one derivative. In most "flat" region starts only 

at fi 2 ~ 10 GeV 2 , where it is hard to believe that cut-off effects are negligible. 
Only for V2b (may be also for O r - 2b ) the window might extend to lower fi 2 . 
It should also be noted that for larger values of /i 2 the scale dependence of 
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Fig. 13. Z RGl for V2 b at (3 = 6.2. 
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Fig. 14. Z RGI for a 2 , at /? = 6.2. 

Z itself (and of i?) is already rather weak so that approximate independence 
of fi 2 is relatively easy to get. Let us now compare operators belonging to 
the same 0(4) representation in the continuum but transforming differently 
under H(4) on the lattice so that their lattice Z's may be different although 
their continuum Z's must coincide. Remember that we have two pairs of that 
kind: V2a , V2b and O r2a , O r2b . The results for V2a lie consistently above 
those for V2b , whereas the separation between O r2a and O r2b is less clear. 
Having an anticipatory look at Table 5 below we note that one-loop lattice 



27 



1.40 



1.35 
1.30 
1.25 



O 
21.20 
N 



1.15 
1.10 
1.05 
1.00 



10" 



o § 



5 °L 



CP 



5 10' 

[i [GeV 2 ] 



Fig. 15. Z RGI for O r2b at /? = 6.2. 
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Fig. 16. Z RG i for O a2 at (3 = 6.2. 

perturbation theory predicts the observed ordering in the case of V2 , 0„ 2 6 
although it underestimates the size of the difference. For O r2a , O r2b , on the 
other hand, it predicts a splitting which is roughly one order of magnitude 
smaller than in the previous case. 

In Figs. 16-18 we plot Zrgi versus the renormalisation scale fi 2 for three op- 
erators containing two derivatives. The size of the numbers as well as their 
spread is larger than in the case of the operators with only one derivative. 
Again, a flat region is observed only for rather large values of /i 2 . One should 
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Fig. 17. Z RGI for O r3 at (3 = 6.2. 
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Fig. 18. Z RGl for a 3 at (3 = 6.2. 

also keep in mind that potential mixing problems in the case of O r3 and V3 
have been neglected. 

In the case of V4 , an operator with three derivatives, it is simply impossible to 
identify a scaling window as Fig. 19 shows. This could be due to the neglected 
mixing problems (cf. Section 2). If this is true, the scaling behaviour should 
be a sensitive test of any attempt to take the mixing into account. But one 
would expect anyhow that cut-off effects might be relatively strong because 
of the large extent of the lattice operator. 
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Fig. 19. Z KGl for O m at /? = 6.2. 



7.^ Comparing (3 = 6.0 and (3 = 6.2 



Let us now compare our results at (3 — 6.0 and (3 = 6.2 for a few representa- 
tive operators. According to the perturbative renormalisation group, the ratio 
of the 2Ts is given by the /i 2 independent factor (5.25) in terms of the bare 
coupling constants g and g' . In Figs. 20-23 we plot Zrgi for (3 = 6.0 and for 
(3 = 6.2 versus the renormalisation scale fi 2 . The corresponding plots for the 
local vector and axial vector currents have already been shown in Figs. 8 and 
9. The f3 = 6.2 data have been rescaled by multiplication with (5.25), which 
indeed moves them closer to the (3 = 6.0 results. Even outside the scaling win- 
dows, the (3 dependence of our data seems to be compatible with the perturba- 
tive expectation (5.25), at least for most of our operators. At first sight, this is 
somewhat puzzling: The perturbative renormalisation group seems to describe 
the dependence on the bare coupling better than the dependence on the mo- 
mentum scale fi. However, the ratio a 2 ((3 = 6.0)/a 2 (/3 = 6.2) « 1.84 = 1.36 2 
is not terribly large and the factor (5.25) differs from one by at most 2 - 3 %. 
Therefore this test of the (3 dependence is not too stringent. 

Still, it is remarkable that the (3 dependence is close to perturbative even for 
rather small values of /i 2 . This may be taken as an indication that the observed 
/i 2 dependence is a real physical effect even in regions where it does not follow 
the perturbative renormalisation group. 
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Fig. 20. Zrqi for O s . The open circles (filled squares) represent the data for [3 = 6.2 
(0 = 6.0). The (3 = 6.2 data have been rescaled perturbatively to (3 = 6.0 by 
multiplication with (5.25). 
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Fig. 21. The same as Fig. 20, but for V2 i 
7. 5 Compilation of results 



In Tables 2 - 4 we present our results for Zq in the MOM scheme at selected 
values of the renormalisation scale. It is to be expected that the systematic 
uncertainties (in particular cut-off effects) are considerably larger than the 
quoted statistical errors. A reasonable measure for the size of the cut-off effects 
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Fig. 22. The same as Fig. 20, but for O vs . 
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Fig. 23. The same as Fig. 20, but for O, 



a 2 - 



might be the width of the band of data points which is typically of the order 
of a few percent for not too small values of /i 2 . 

In Table 5 we give the nonperturbative values for Z in the MS scheme to- 
gether with Zq and for /i 2 = a -2 . In order to determine the nonpertur- 
bative Zo& for this value of fi 2 we have interpolated linearly in In /i 2 between 
the two neighbouring results. The given error is the larger of the two statis- 
tical errors. Disregarding the results at other values of /z 2 we see that, with 
the exception of the operators containing one covariant derivative, tadpole 
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Table 2 

Renormalisation factors in the MOM scheme at selected values of the renormalisa- 
tion scale /j, 2 a 2 in lattice units. The errors are purely statistical. 





fi 2 a 2 


qA 




O p 






0.318 


0.809(2) 


0.645(2) 


0.229(3) 


0.602(9) 




0.626 


0.779(2) 


0.6658(6) 


0.360(3) 


0.658(4) 




0.935 


0.7803(9) 


0.6814(8) 


0.431(2) 


0.679(1) 




1.320 


0.7816(7) 


0.6893(7) 


0.487(1) 


0.720(1) 


6.0 


1.860 


0.7850(4) 


0.7024(6) 


0.5347(7) 


0.742(1) 




2.477 


0.7878(4) 


0.7088(5) 


0.5772(5) 


0.7752(4) 




3.942 


0.7995(2) 


0.7286(3) 


0.6320(2) 


0.8041(2) 




5.330 


0.8094(2) 


0.7543(3) 


0.6731(2) 


0.8100(2) 




7.797 


0.8268(2) 


0.7761(3) 


0.7127(2) 


0.8352(2) 




0.313 


0.7954(8) 


0.6951(3) 


0.357(4) 


0.605(7) 




0.587 


0.797(1) 


0.7075(3) 


0.463(3) 


0.6772(9) 




0.930 


0.7976(4) 


0.7189(3) 


0.524(1) 


0.709(3) 




1.272 


0.7994(4) 


0.7251(3) 


0.566(1) 


0.738(2) 


6.2 


1.855 


0.8022(4) 


0.7352(2) 


0.6078(9) 


0.7646(8) 




2.369 


0.8053(2) 


0.7435(1) 


0.6328(4) 


0.7790(7) 




4.014 


0.8161(1) 


0.7649(1) 


0.6877(3) 


0.8100(4) 




5.385 


0.8254(1) 


0.77760(9) 


0.7150(1) 


0.8289(2) 




7.921 


0.84238(8) 


0.80054(5) 


0.75556(5) 


0.8542(2) 



improvement works in the right direction. In the case of the local vector and 
axial vector currents it even overshoots somewhat. However, only in these 
latter cases the improvement is quantitatively satisfactory. 

The renormalisation factors for the operators without derivatives (O s , O p , 
O v , 0£) have been calculated with similar nonperturbative methods by other 
groups, too. However, most of these studies use either an improved fermionic 
action or work at different values of /3 so that a direct comparison is impossible. 
We are aware of only one other paper dealing with Wilson fermions at /3 — 6.0 
and (3 = 6.2 [28]. Let us compare our results as given in Table 2 with Table 3 
of Ref. [28]. 
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Table 3 

Renormalisation factors in the MOM scheme at selected values of the renormalisa- 
tion scale /j, 2 a 2 in lattice units. The errors are purely statistical. 









^2,6 


V3 


o V4 




0.318 


1.467(8) 


1.430(9) 


2.42(3) 


5.7(1.3) 




0.626 


1.352(9) 


1.294(9) 


2.08(1) 


2.94(2) 




0.935 


1.219(3) 


1.160(3) 


1.692(7) 


2.44(1) 




1.320 


1.162(2) 


1.139(3) 


1.596(3) 


2.31(2) 


6.0 


1.860 


1.1042(9) 


1.069(1) 


1.454(2) 


2.000(4) 




2.477 


1.0999(7) 


1.0486(5) 


1.431(2) 


1.879(8) 




3.942 


1.0366(4) 


1.0527(8) 


1.3167(9) 


1.701(3) 




5.330 


1.0254(2) 


1.0094(2) 


1.2753(5) 


1.64(1) 




7.797 


1.0271(2) 


1.0107(2) 


1.2615(4) 


0.909(2) 




0.313 


1.408(5) 


1.304(2) 


1.98(2) 


3.02(5) 




0.587 


1.236(3) 


1.210(7) 


1.731(8) 


2.53(6) 




0.930 


1.158(2) 


1.128(4) 


1.529(2) 


2.084(5) 




1.272 


1.128(1) 


1.102(3) 


1.4696(6) 


1.950(3) 


6.2 


1.855 


1.0774(7) 


1.057(2) 


1.3709(5) 


1.799(3) 




2.369 


1.0583(5) 


1.038(1) 


1.3319(4) 


1.73(1) 




4.014 


1.0313(2) 


1.0190(7) 


1.2680(4) 


1.583(7) 




5.385 


1.0321(2) 


1.0181(5) 


1.2543(3) 


1.466(2) 




7.921 


1.0210(2) 


1.0211(4) 


1.2296(5) 


1.463(1) 



• In the case of the local vector and axial vector currents we find some discrep- 
ancies, which, however, decrease with increasing fi 2 . At (3 = 6.2 we obtain 
consistent results for the larger values of fi 2 . The deviations could be due to 
the fact that we use only the transverse components of the currents. Indeed, 
averaging over the four Z's determined from the four components of the lo- 
cal vector current by the method of Section 3, we arrive at Fig. 24, where 
also the results of Ref. [28] are shown. We observe satisfactory agreement 
as well as a relatively large spread of the data, which is not present in our 
method based on the transverse components (see Fig. 8). 

• The pseudoscalar density shows a different behaviour. Our Z's are consis- 
tently smaller for (3 = 6.0, but always larger for (3 = 6.2 (although the 
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Table 4 

Renormalisation factors in the MOM scheme at selected values of the renormalisa- 
tion scale /j, 2 a 2 in lattice units. The errors are purely statistical. 





O a2 


o n , a 


Or Xb 


O r3 


0.318 


2.44(2) 


1.539(5) 


1.60(1) 


2.52(3) 


0.626 


2.06(1) 


1.382(7) 


1.417(10) 


2.19(1) 


0.935 


1.796(7) 


1.259(4) 


1.235(3) 


1.724(7) 


1.320 


1.636(5) 


1.189(3) 


1.211(3) 


1.629(4) 


6.0 1.860 


1.503(3) 


1.129(1) 


1.1230(9) 


1.478(2) 


2.477 


1.494(3) 


1.1162(7) 


1.1040(5) 


1.455(2) 


3.942 


1.360(1) 


1.0458(4) 


1.0942(8) 


1.331(1) 


5.330 


1.2852(6) 


1.0282(2) 


1.0387(3) 


1.2869(6) 


7.797 


1.2773(6) 


1.0308(2) 


1.0363(2) 


1.2710(5) 


0.313 


2.17(1) 


1.457(6) 


1.380(3) 


2.01(1) 


0.587 


1.776(5) 


1.265(3) 


1.281(8) 


1.760(8) 


0.930 


1.594(2) 


1.181(2) 


1.183(4) 


1.548(2) 


1.272 


1.5152(8) 


1.145(1) 


1.154(3) 


1.4884(8) 


6.2 1.855 


1.4124(5) 


1.0938(9) 


1.098(1) 


1.3869(5) 


2.369 


1.3564(3) 


1.0704(6) 


1.0737(10) 


1.3455(4) 


4.014 


1.2782(3) 


1.0358(3) 


1.0460(6) 


1.2780(3) 


5.385 


1.2479(3) 


1.0351(2) 


1.0417(4) 


1.2627(3) 


7.921 


1.2241(3) 


1.0194(2) 


1.0386(3) 


1.2372(4) 



differences are not large). Here the chiral extrapolation might be responsi- 
ble for the discrepancies. Especially for the lower values of /i 2 a 2 we observe 
deviations from a linear 1/k dependence, which make the extrapolation 
somewhat questionable and ambiguous. 
• For the scalar density we find satisfactory agreement if /i 2 a 2 is not too small. 
Like in our data, the ratio Z p/Z s does not become constant even for the 
largest values of /i 2 considered in [28]. 

In all cases one should keep in mind that the lattice used in [28] at /3 — 6.2 
is smaller than ours. Furthermore, the momenta may differ in direction even 
if they are close in length, hence part of the differences could also be lattice 
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Table 5 

Renormalisation factors in the MS scheme for (3 = 6.0 at fj? = a~ 2 = 3.8 GeV 2 and 
for (3 = 6.2 at /U 2 = cT 2 = 7 GeV 2 .The errors of the nonperturbative values Zq are 
purely statistical. 



o 






ypert 
O 






Zq 


ypert 
O 




3 


= 6.0 




3 




6.2 






1 1332(291 


0.9731 


9892 


1 0732(19") 


0. 


.9759 


9895 


O 


1 0836(24) 


0.9462 


0.9784 


1 0426(36) 


o. 


.9518 


0.9791 


O 


1 4973(67) 


1.1933 


1.1024 


1 3458(15) 


1 


.1778 


1.0991 


O 


2 0042(121) 


1.5021 


1.2299 


1 7497(85) 


1 


.4565 


1.2225 


(0 


1 54D9(64) 


1 193D 


1 1D93 


1 3969(16) 


i 

± . 


1 776 

ill IU 


1 099D 


Or %a 


1.1615(33) 


0.9927 


0.9971 


1.0956(18) 


0. 


.9935 


0.9972 


Or 2 , b 


1.1541(24) 


0.9965 


0.9986 


1.0912(36) 


0. 


.9968 


0.9986 


O r , 


1.5292(68) 


1.2108 


1.1086 


1.3628(15) 


1 


,1934 


1.1051 


Q S 


0.7718(16) 


0.8209 


0.8906 


0.7897(29) 


0. 


.8337 


0.8942 


o p 


0.4934(24) 


0.6430 


0.8092 


0.5888(15) 


0. 


.6731 


0.8154 


qV 

M 


0.6833(8) 


0.6795 


0.8259 


0.7206(3) 


0. 


.7060 


0.8315 


qA 


0.7821(9) 


0.7684 


0.8666 


0.7978(4) 


0. 


.7864 


0.8709 



artifacts. 

The renormalisation constant of the local vector current can also be computed 
from its correlation functions with hadron sources by comparing with the 
corresponding correlation functions of the conserved vector current, whose 
renormalisation constant is known to be one. A collection of results obtained 
with various versions of this method at f3 — 6.0 is given by Sachrajda [29] (cf. 
also [30]), albeit without extrapolation to the chiral limit. Since the quark mass 
dependence is rather mild, we may nevertheless compare with our numbers as 
displayed in Fig. 8. The values given in Ref. [29] cluster around 0.73 (computed 
from three-point functions) and 0.57 (computed from two-point functions). 
Sachrajda attributes this large difference to discretization effects. This is in 
accordance with our interpretation of the strong scale dependence of the Z 
for the local vector current (see Subsection 7.1). More recently, the JLQCD 
collaboration has applied this method with Wilson fermions at f3 — 5.9, 6.1, 
and 6.3 [31]. Interpolating their results [32] one finds 0.57 at (5 — 6.0 and 0.62 
at = 6.2. 
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Fig. 24. Z for the local vector current at (3 = 6.0. The filled squares represent our 
data averaged over the four components of the current. The open squares are results 
taken from Ref. [28]. 

The renormalisation of the local vector current and the axial vector current 
can also be studied by means of the Ward identity method [33] . However, the 
values obtained for Wilson fermions at (3 = 6.0 by Maiani and Martinelli [30] 
(0.79(4) for the vector current and 0.85(7) for the axial vector current) refer to 
point-split lattice currents and are therefore not immediately comparable with 
our numbers. On the other hand, for the ratio Z v/Z a of the renormalisation 
constants of the local currents they give the values 0.85(5) and 0.79(5) which 
are lower than (though not completely incompatible with) our result of 0.87 
obtained at /i 2 a 2 = 1. From the Z's which the JLQCD collaboration has 
calculated for the axial vector current by the Ward identity method at (3 = 5.9, 
6.1, and 6.3 one gets by interpolation [32] approximately 0.75 (f3 = 6.0) and 
0.77 (f3 = 6.2), which is somewhat lower than our values for /i 2 a 2 = I. 

There have been various proposals to increase the accuracy in the determi- 
nation of renormalisation constants, in particular by using modified versions 
of tadpole improvement. For Z v, Z a, and the ratio Z p/Z s these are 
discussed in Ref. [34], where also a comparison with results from the Ward 
identity method is given. 



8 Conclusions 



We have presented a comprehensive study of nonperturbative renormalisa- 
tion for various operators in the framework of lattice QCD. In particular, 
twist-2 operators appearing in unpolarised (polarised) deep-inelastic scatter- 
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ing have been studied for all spins < 4 (< 3). We worked with standard Wil- 
son fermions at two different lattice spacings. Due to the use of momentum 
sources we achieved a high statistical accuracy so that systematic uncertain- 
ties (cut-off effects, in particular) are clearly visible. In order to make contact 
with perturbation theory one would like to identify scaling windows where the 
renormalisation scale fi is large enough to make perturbation theory trustwor- 
thy but small enough to allow to neglect cut-off effects. In this respect some 
operators did not follow our (naive) expectations. Whereas the axial vector 
current and the scalar density behave as expected, the local vector current and 
the pseudoscalar density do not show perturbative scaling. A second thought, 
however, revealed that there might be good reasons for this behaviour re- 
lated to the spontaneous breakdown of chiral symmetry and there is nothing 
fundamentally wrong. 

On the other hand, the operators containing one or more covariant deriva- 
tives seem to approach perturbative scaling only for rather large values of the 
renormalisation scale where cut-off effects might influence the results. Can we 
understand these deviations from perturbative scaling which seem to persist 
up to fj, 2 ~ I OGeV 2 ? If the observed scaling with (3 is not only an artifact of our 
two P values being relatively close to each other, the deviations from scaling 
in n should be interpreted as real physics, though in a finite volume. (Recall 
that our lattices at the two f3 values have roughly the same physical size.) So 
one possibility would be finite-size effects, although the experience with other 
observables suggests that they should be small. But genuinely nonperturba- 
tive effects cannot be ruled out, and one might also think of renormalons. 
In any case, one does not feel very comfortable when combining these non- 
perturbative Z's with perturbative Wilson coefficients using a scale of, say, 
^ = 4GeV 2 . 

One remedy would be a nonperturbative calculation of the Wilson coefficients. 
But one can also try to reach larger values of fi 2 . In order to avoid strong 
cut-off effects one would then have to work at smaller lattice spacings, i.e. 
at /3 > 6.2. However, this raises a new problem. At present, we have been 
careful to work with lattice sizes L large enough that we believe finite-size 
effects are not serious. But if we try to reach much smaller a's it soon becomes 
impossibly expensive to keep the physical lattice size sufficiently large. So we 
need a reliable method of measuring Z nonperturbatively, even when finite- 
size effects are not negligible. Such a scheme has been suggested in Ref. [35]. 
The key point is to note that although the Green functions depend on L, the 
renormalisation constants do not. This means that we can nonperturbatively 
measure the ratios of Z's at different lattice spacings by looking at the ratios 
of the corresponding bare Green functions, both measured on lattices of the 
same physical size, and with bare masses chosen so that the renormalised 
mass is the same in both cases. Note that to apply this scheme, we need good 
knowledge of the lattice spacing as a function of f3, which may be a problem. 
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Further improvement could come from continuum perturbation theory. In the 
case of the scalar density the use of the two-loop formula for Z^om an d of the 
three-loop anomalous dimension improves the (approximate) scale indepen- 
dence of -Zrgi- For the operators with derivatives we only have the two-loop 
anomalous dimension and a one-loop calculation of ^mom- Recently, however, 
first steps towards a three-loop computation of the anomalous dimension have 
been undertaken [36]. From the results already obtained, a two-loop expres- 
sion for Z^om can be extracted, unfortunately only in the Feynman gauge, 
whereas we would need the Landau gauge. 
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